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Abstract 

Let K be a field of cliaracteristic zero, a algebraic of degree n over K. 
Given a proper parametrization i/) of a rational curve C, we present a new 
algorithm to compute the hypercircle associated to the parametrization 
ip. As a consequence, we can decide if C is defined over K and, if not, to 
compute the minimum field of definition of C containing K. The algorithm 
exploits the conjugate curves of C but avoids computation in the normal 
closure of K(a) over K. 

1 Introduction 

Let K(a) be a computable characteristic zero field with factorization such that 
K is finitely generated over Q as a field and a is of degree n over K. 

Let ipi^) = (^"1(01 ■ • ■ I V'm(^)) be a proper parametrization of a rational 
spatial curve C, where ipi G K(a)(t), 1 < i < m. The reparametrization problem 
ask for methods to decide in C is defined or parametrizable over IK and, if 
possible, compute a parametrization of C over K. 

In [T] , the authors proposed a construction to solve this problem introducing 
a family of curves called hypercircles and avoiding any implicitization technique. 
Starting from the parametrization ip, they construct an analog to Weil descente 
variety to compute a curve U called the witness variety or the parametric variety 
of Weil. This curve exists if and only if C is defined over K and we can obtain 
a parametrization of C with coefficients in K easily from a parametrization of U 
with coefficients in K. Efficient algorithms to compute a parametrization of U 
with coefficients in K are studied in provided we are able to find a point in 
U with coefficients over K. 

The definition of U is done under a parametric version of Weil's descente 
method. In the proper parametrization ip — (^i: • • ■ ,4'm), ipi G JK(Q!)(i) with 
coefficients in IfC(a), we substitute t = J27=o Q^^^ij where to, ■ ■ ■ ,in-i (where n 
is the degree of a over K). 

We can rewrite: 




- a^Xij (to, • • . , ^n-i), Ay — -jj- e K(io, • • ■ , in-i) 

1=0 

In this context we have the following definition: 
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Definition 1. The parametric variety of Weil Z of the parametrization ip is he 
Zariski closure of 

{F,j ^0 \ 1 <i <n-l, I < j < N}\{D ^0} C F". 

Much is known about Z, it is always a set of dimension or 1. It is of 
dimension one exactly in the case that C is defined over K (See [I], [1]). In 
this case, Z contains exactly one component of dimension 1 that is the searched 
curve U. 

The computation of the curve U from its definition is unfeasible except for 
toy examples. The curve U is defined as the unique one dimensional com- 
ponent of a the difference of two varieties A — B. This already is a hard 
enough problem to look for alternatives, but this method also uses huge poly- 
nomials. If i'iit) = ni{t)/d{t) and d = d{a,t) G K[a,t]. Let M{x) be the 
minimal polynomial of a over K. In the generic case, the denominator D is 
D = i?es^(d(z, ^"Jq"^ M(z)) which is typically a dense polynomial of de- 
gree dn in n variables. Hence, the number of terms of the polynomial D alone 
is not polynomially bounded in n. 

The aim of the article is to present an algorithm to compute the variety U 
that is polynomial in d and n and, if C is not defined over K, to compute the 
smallest field L, K C L C K{a) that defines C. The article is structured as 
follows. First we introduce in Section [5] the geometric construction that will 
allow us to derive an efficient algorithm. Then, we show in Section |3] how to 
compute efficiently some steps of the algorithm. Last, in Section SI we study 
the complexity of the algorithm and some running times comparing with other 
approaches. 



2 Synthetic construction of Hypercircles 

The problem of parametrizing C over K can be translated to the problem of 
parametrizing U. In the case that C can be parametrized over IK, then U is a 
very special curve called hypercircle. 

Definition 2. Let G K{a){t) represent an isomorphism of¥{t), a,b,c,d G 
K(a), ad -6c 7^0. Vriie 

= Xo{t) + aAi (i) + • ■ • + a"-' A„_i (t) 

ct + d 

where Xi{t) £ K.{t). The hypercircle associated to for the extension IK C 

K(a) is the parametric curve in F" given by the parametrization (Aq, . . . , A„_i). 

If C cannot be parametrized over K and K is small enough (that means that 
it is finitely generated over Q as a field, that we can always assume without loss 
of generality), then there always exists an element /3 algebraic of degree 2 over IK 
such that [IK(/3, a) : K{a)] = n and C can be parametrized over K(/3), see [12] for 
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the details. In this situation is a hypercircle for the extension ]K(/3) C K(q, /?). 
That is, there is an associated unit f|^, but with a, b,c,d& K(/3). 

Thus, the curve U is always a hypercircle for certain algebraic extension. 
So all the geometric properties of hypercircles studied in [B] hold for U except, 
maybe, the existence of a point in U H K". We will exploit the geometric 
properties of hypercircles to derive our algorithm. We start with the fact that 
U is always a rational normal curve in F" defined over K (See [6]) and the 
synthetic construction of rational normal curves as presented in [5]. 

Let us recall the construction of conies by a pair of pencil of lines. Let £,{t) 
and ^{t) be two different pencils of lines in the plane with two different base 
points /q ^ /o and let C be a conic passing trough Iq and /q. Then, C induces 
an isomorphism u : £,{t) ^{t) given by extending the map w(£(to)) — 5^(so) if 

£{to)nc~{io}^d{so)nc-{fo}. 

Conversely, an isomorphism u between £,{t) and ^{t) defines a line or a conic 
passing through the base points. There is a proper parametrization of this curve 
given by iK^£(i)n;?(u(i)). 

Example 3. Let C = + — 1 be the unit circle. And take the pencils of 
lines that passes through the points at infinity of the circle [1 : i : 0], [1 : — i : 0] . 
£(t) ^{x + iy^ t}, g-(i) = {x-iy = t}. In this case C n Z{t) = ^^^) 

and C r\^{t) = (^^^, ^*'2t^ )- ^'^ case, the isomorphism between the pencils 
is given by u(t) ~ 1/t. Now, let us take the isomorphism u{t) ~ {t + i)/t. Then, 
the conic defined by u from the two pencils of lines is x'^ + y^ — x — iy — i. Which 
is a conic passing through the base points, although not defined over Q. 

More generally, the same geometric construction applies to rational nor- 
mal curves of degree n > 2 in F" as explained in We only show the 
special case of this construction that is relevant for hypercircles. U U is a 
hypercircle, it is known that 14 can be parametrized by the pencil of hyper- 
planes £o = {J2i=o ~ i^.- '^^i^ pencil of hyperplanes yield to a proper 
parametrization </> = ((/)o(t), . . . , (j)n^i(t)) of the hypercircle with coefficients in 
K(a) that is called the standard parametrization of the hypercircle and has been 
studied with detail in 7 . Since the hypercircle is always a curve defined over 
K, it is invariant under conjugation and it can also be parametrized by the 
conjugate pencil of hyperplanes. 

Let us fix some notation. Let a — ao, ai, . . . , a„_.i, the conjugates of a over 
K in F. Let ai,0<i<n — Ihe K- automorphisms of F such that ai{a) = ai 
and (To — Id. If we have a rational function f{t) e K(a)(xi, . . . , Xr), we denote 
by f^^ — Cjif) £ K{aj){xi, . . . ,Xr) that results applying Cj to the coefficients 
of /. If C is the original curve, then we denote by the conjugate curve 
C"^ = {a{x)\x £ C}, where c^x) is applied component-wise. is clearly a 
rational curve with proper parametrization . 

It is known 2] that C is defined over K if and only if C^C^' l<i<n — 1 
if and only if tjj'^^ parametrizes C, I < i < n — 1. 
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The conjugate pencil of hyperplanes £>j{t) ~ = 1 ^ J n — 1 

also parametrizes U, yielding the conjugate parametrization (j)"' (t) ~ Cj 

The hypercircle then induces an isomorphism Uj(i) between £o(0 ^ji^) 
given by (£o(io)nZ//) — -ff = {2j{uj{to))r\U)—H for all but finitely many param- 
eters to, where H is the hyperplane at infinity of P(F)". So = (j)'^' {uj{t)), 
from which Uj{t) — {(f>'^^)~^ o 0. But, by construction, (^'^^ )~^(a;o, . . . ,Xn~i) = 
127=0 ^i^d Uj{t) = X]"=o '^)'Piit)- Conversely, a set of isomorphisms uj : 
2o{t) ^ ^j{t), < j < n — 1, uo{t) = t, defines a rational normal curve given 
by ^ ~^ n"=o ^6 '^^^ recover the standard parametrization of the 

hypercircle if we know the isomorphisms Uj, < j < n~l, where uo{t) — t. The 
standard parametrization </> is the unique solution of the Vandermonde linear 
system of equations: 

/ 



a2 



/-I 



\4>n-l) 



/uo{t) = t\ 
Ul{t) 



(1) 



with coefficients on the normal closure of a over K. 

As in the planar case, if the automorphisms are generic enough, the curve U 
will be of degree n. In this case we say that W is a primitive hypercircle. There 
may be cases in which the curve U is of degree less than n. If this is the case, 
the degree of U must be a divisor of n and is related with the field of definition 
of the place of C corresponding to 4'{t = oo), as showed in |12) . 

The good news is that we can compute easily the automorphisms Uj (t) from 
the parametrization 'ijj{t) alone. 

Theorem 4. Let ^p{t) G ]K(a)(i)'" be a proper parametrization of C and as- 
sume that C is defined over K. Let (j)(t) be the standard parametrization of the 
associated hypercircle lA. Let Ui be a automorphism of¥. Let 0'^' — CFi{(j)), 
ip"'' = (Ti(V') be the conjugate parametrizations and u„.{t) — {(jf^)^^ o (j) be the 
conjugation isomorphism induced by U in the pencil of hyperplanes £o and £0-; ■ 
Then Ua-^ = o -0. 

Proof. We identify C with the diagonal curve A in the variety C x C^^ x . . . x 
C°'""i, A = {{x, . . . ,x)\x e C} the hypercircle W is a curve such that the map 



U C X X ... X C""-! 

(n—1 n—1 n—1 \ 

j=0 j=0 j=0 J 

Is a birational map between U and A = C. See '2 for the details. This 

means that V'lX^j^o ^j'^"') ~ '^'^'{J2]=o ^j'^l) for the points of the hypercircle. 
If we plug the standard parametrization of the hypercircle in this equality, we 
get that 

n—1 n—1 
j=0 j=0 
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From which u^^ 



-1 



o ip. 



□ 



Hence the isomorphism Uj induced by the hypercircle in the pencil of hy- 
perplanes £{t) and is the change of variables needed to transform the 

conjugate parametrization 'ip'^^{t) into ip{t). We can compute Uj using gcd. 

Theorem 5. Letijj^{t) = n,{t)/d,{t) andipyit) = nl' [t) / (f^' {t) be the numera- 
tors and denominators ofijj andil^K Then, if C is defined over K, the numerator 
of s — Uj (t) is a polynomial of degree 1 in t and in s that is the common factor 
of the set of polynomials 



On the other hand, if C is not defined over K, there is an index 1 < j < n ~ 1 
such that gcd{Ba-j) — 1. 

Proof. This result follows directly from the geometric interpretation. First, 
assume that C is defined over K. It is clear that the numerator of s — Uj (t) is a 
common factor of the set Bc^ . Let f{t, s) be the gcd of and let p — tp{to) G C 
where to is a generic evaluation of t. The roots of f{to,s) are solutions of the 
system of equations '0^ ' (s) — Pi- But, since ip'^^ is birational, for all but finitely 
many to there is only one solution, {ip'^^ (p) . Hence, the degree of / with 
respect to s is one. By symmetry, the degree of / with respect to t is also one. 
It follows that / must be the numerator of s — u{t). 

Now, assume that C is not defined over K. Then, there is an index j such that 
C ^ . In this situation, for all but finitely many evaluations t — to, the system 
of equations ip'''{s) — ^{to) has no solution. It follows that gcA{B„.) = 1. □ 

So, we can compute K-definability and the standard parametrization of the 
hypercircle lA by the following method: 

• For each conjugate aj, Compute a{t) + sh{t), the gcd of ni{t) ■ d^' {s) — 

(s) ■ di{t),l < i < m. If one of the gcd is one, then the curve is not 
defined over K and we are done. 

• Set Uj = -a{t)/b{t). 

• Solve the linear system of equations ([T]) whose coefficients are rational 
functions in t with coefficients in the normal closure of K(a). 

However, computing these bivariate gcd are expensive and, moreover, in the 
worst case, we will have to solve a linear set of equations with coefficients in an 
extension of K of degree nl. Next section address the problem of how to perform 
this algorithm efficiently. 

3 Efficient Computation of the Hypercircle 

We have shown how to compute (t) by computing the gcd of the polynomials 
in Ba- We already now that, if C is K-definable, the gcd has degree 1 in t and 



B,^ = {n,{t) ■ d^^ is) 



n'^'is) ■d^{t),l <i<m}. 
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s, so the best suited algorithms for computing the gcd seem to be interpolation 
algorithms. Since we are only interested in Ua- and this linear fraction is an 
automorphism of P^(F), we only need to know the image of three points to,ti, t2 
under u^. From Theorem [SJ for almost all ti, Ua{ti) = (si) if and only if 
ip{ti) = ip°'{si). Hence, each Si is the common root of the polynomials: 

■dj(s) -nj(s),l <i<m 

Si that can be computed by means of gcd of univariate polynomials in 
K{a,a{a)). 

If C is defined over K then only finitely many parameters tk will fail to 
provide a valid Sk- Essentially the parameters tk can fail if ip{tk) is a singular 
point of the curve or if it cannot be attained by a finite parameter s by the 
parametrization ip"^ . 

On the other hand, if C is not defined over K, then there is an automorphism 
a such that C ^ C" . For this permutation, there are only finitely many parame- 
ters tk such that ijj{tk) G C n . Hence, if we want to follow this approach and 
do not depend on probabilistic algorithms that may fail or give wrong answers, 
we need bounds to detect that the curve is defined over IK or not. 

Theorem 6. Let C C F™ be a rational curve of degree d given by a parametriza- 
tion ip G (IK(a))™. Let ai be any conjugate of a over K. Take ti,...,tk G F 
parameters then: 

• If C is definable over K, then we can compute Ui from three correct solu- 
tions of the system of equations ip'^^s) = tk- 

• If C is defined over K, then at most d'^ — 2d + n + 1 parameters can fail to 
give a correct answer. 

• If C is not defined over K, then at most d^ parameters tk will give a fake 
answer Sk ■ 

Proof. We have to compute the inverse of the point ip{tj) under the parametriza- 
tion ?/''^'(s). For each tk, this computation is done using univariate gcd. If we 
want to restrict to affine points, we have to eliminate d potential parameters of 
the denominator of iJj. Then, for an affine point ipitj), there can only be one 
point that is not attained by a finite parameter of tp"^ . Since we have n — 1 
possible conjugates, then there may be n — 1 points that are not attained by 
a finite parameter in one of the conjugate parmetrizations. So, if we get two 
different parameters tj such that is well defined but that ■0°'*(s) = tj have 

no solution (the corresponding gcd is 1), then the curve is not defined over K. 
Now, it may happen that the gcd is of degree > 1 . This can only happen if the 
point is singular in C. Since C is of genus and degree d, it can have at most 
{d — l){d — 2)/2 singularities. The number of different parameters whose image 
is a singularity is maximal if every singularity is ordinary. We have to maximize 

multp{C) 

p^sing{C) 
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subject to 



multp(C){multp{C) -l)^{d- l){d - 2) 

See [5] Theorem 2.60 for details. But clearly, for any singular point multp{C) < 
mult p{C) {mult p{C) — 1) So X]pGsmg(c) ™"'*p(^) — (^^ ~ — 2) and the 
equality is attained if every singularity is an ordinary double point. 

Thus, the maximal number of parameters that cannot be used to compute 
Ui is bounded by d parameters corresponding to the points at infinity plus 
n — 1 points that might not be attained by a finite parameter in a conjugate 
parametrization ij:'^ plus {d — \){d — 2) parameters whose image are singular 
points. This gives the bound d^ — 2d + n + 1. 

Suppose now that C is not defined over K. Let at be such that C ^ C^'. A 
parameter gives a fake answer for computing Ui if "(/"(^o) is smooth in C"^' and 
is attained by a unique parameter sq by V''^'- But, by Bezout, CCiC^ contains at 
most different points. So, there can be at most d^ such bad parameters. □ 

Remark 7. In order to check that a parameter tk is a good parameter or not 
we can do the following: 

• Iftk is a root of the denominator of tp, then tk is a bad parameter. 

• //gcd(?/'(ifc) • dfls) — n"{s), 1 < i < m) = 1, then it is a bad parameter. It 
is a point that is not attained by the parametrization -0°" . If C is defined 
over ]K there can be at most one bad parameter that happens to be in this 
case that corresponds to ijj'^it = oo). 

• If deg{gcd{ip{tk)-df {s)—nf {s),l <i < m)) > 1 thentk is a bad parameter, 
since 'ip{tk) is a .singular point. 

• //deg(gcd(V'(ife) -dHs) -nf(s),f < j < m)) = f but ip{tk) ^ ip" (oo) then 
tk is a bad parameter, ^p(tk) is .singular. 

• // deg(gcd(V'(ife) • df{s) - nf{s), 1 < i < m)) ^ 1 and V(ife) / ^'''(oo) 
then tk is a good parameter, we compute Sk solving the linear equation in 
s given by the gcd. 

Hence, if C is defined over IK, we can compute each Uj by interpolation. We 
will need at most d^ — 2d + n + l + 3 = d'^ — 2d + n + 4: parameters. In practice 
however, we will almost always need only 3 parameters. Note also that if we 
choose the parameters in K, then all computations needed to compute Ui are 
done in K(a, ai), that is a extension of degree bounded by n(n — 1). 

If C is not defined over K it can happen two things while trying to compute 
Uj. With high probability, we may find two different parameters such that 
gcd(V'(ifc) • df{s) — nf{s),l < i < m) — 1 and this is a certificate that the 
curve is not defined over K. On the other hand, we may succeed computing 
Uj . This may happen if C = for this specific aj or if we have chosen three 
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parameters tk such that tp{tk) E C D C"'^ . So, if we have computed all the linear 
fractions Uj{t) but we want a certificate that C is defined over K, we only need 
to check that tpit) = 4>'^' {uj{t)), 1 < i < n. In the case that computing this 
composition may be expensive, we can try to check the equality evaluating in 
several parameters t. tp{t) and ip'^^{uj{t)) are rational functions of degree d, so 
if they agree on 2d + 1 parameters where both parametrizations are defined, 
then ■ijj{t) = ■ijj'^^ {Vjit)) and C = . But there are d parameters where C is not 
defined and other d where is not defined. So, if we want a certificate that 
C = by evaluation, we will need to try at most 4d + 1 parameters in the 
worst case. So (n — l)(4rf+ 1) evaluations to check all conjugates. 

Example 8. Let us show that the bounds given can be easily proven to be sharp 
if we allow the parameters to be in F and d>n. Let K(a) be normal over K of 
degree n and ai, . . . , cr„_i be K- automorphisms that send a onto its conjugates. 
The common denominator of the parametrization of the curve will be g = (t + 
1) . . .{t + d) so that the parameters — 1, . . . , —d will fail in the algorithm. Let 
us write a component as f{t) = {af^ + ad-it'^~^ + . . . + a\t + ao)/g{t), where 
the cii are indeterminates. Impose the conditions f{i) = o"i(a). This is a linear 
system of equations in the ai representing an interpolation problem. We have 
n — 1 conditions and d unknowns in the system and n — l<d. Hence, there are 
infinitely many solutions to the system and we can take two generic solutions 
fi{t),f2{t). The curve tp{t) = {fi{t), f2{t)) will fail to give a correct answer for 
t = —1, . . . , —d due to the denominator and for t = i, i = 1, . . . ,n — 1 because 
ijj(t — oo) = (a, a), so tp'^it = oo) = {a{a.), ^{a)) = ip{i). Finally, if we have 
chosen /i,/2 generic, the only singularities of ^ will be simple nodes in the 
affine plane. Thus, there will be {d — l)(c? — 2) parameters that will yield to a 
singularity. 

For a specific example, take K = Q, a a primitive 5 — th root of unity, so 
that n = 4. Let the degree be d = A. If we perform the construction above, we 
get the relations in the coefficients of f: 

ao = -603 + (1440a^ + lOSOa^ + IQAAa + 1920) 

ai = lias + (-17400!^ - lAAQa^ - 1380a - 2700) 

02 = -6a3 + (420a^ + 360a^ + 335a + 780) 

If we compute fo and fi substituting 03 by and 1 respectively, we get the 
parametrization 4>{fo,fi) of a rational curve of degree 4 with three nodes, such 
that the nodes are attained by the roots of t^ + (-420a^ + 60a^ - 90a - 102)t^ + 
(-59220a3-171720a2+85110a-214952)i4+(688980a3+1237740a2-450750a+ 
1759626)t3 + (-2309580a3 - 3135240a2 + 714450a - 4869077)^^ + (2877600a3 + 
3308280a2 - 391560a + 5387628)t + (-1197360a3 - 1231920a2 + 42840a - 
2063124). 

Now, we show how to avoid in some cases some computations of Uj using 
conjugation. 
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Proposition 9. Assume thatC is defined owerK. Letai ^ Uj he two conjugates 
of a over K. Suppose that ai, aj are also conjugated over K(a) and that t is a 
M^(a) -automorphism o/F such that T{ai) = aj. Then T{ui) = Uj. 

Proof. All operations to compute Uj are evaluating rational functions with co- 
efficients in K(a, ai) at parameters in IK (or even Z), compute gcd of univariate 
polynomials with coefficients also in K(a,ai) and solving a linear system of 
equations. These operations commute with conjugation by a. Thus, if r is a 
K(a)-automorphism such that T{ai) — aj, we can conjugate by r at every step 
of the method to compute Ui. Hence, T{ui) = Uj. □ 



If the Galois group of K(a) over K is the permutation group Sm we will only 
need to compute one automorphism Ui making computations in a number field 
of degree n{n — 1). On the other extreme, if K C IK(a) is normal, we will have 
to compute n — 1 different automorphisms Ui, but the computations will be in 
the smaller field K(q;). 

Now, we show how to avoid computing in the normal closure of K(q;) over 
K to solve the linear system of equations [TJ This system is given by a Vander- 
monde matrix, so we are dealing with an interpolation problem. If the standard 
parametrization searched is ((^o, • ■ • , 4>n~i)- Then, the polynomial 

F{x) = (/)o + + . . . + (/)„_ia;""i e K{a){t)[x\ 

is the unique polynomial of degree at most n — 1 such that F{ai) — Ui{t), 
0<i<n— 1. F can be computed by Lagrange interpolation 

{x - ap) . . .{x - ai-i){x - at+i) . . . {x ~ Qn+i) , , 
^ [ai - ao) . . . [ai - ai-i)[at - aj+ij . . . [ai - an-i) 

Let us take a look at each term: 

{x - ao) ...[x - ai_i)[x - a^+i) . . . [x - a„_i) 



(ai - ao) . .. (ai - q;i_i)(q;j - a^+i) . . . [a^ - a„_i) 

The numerator is M[x)/[x — ai) = m[ai,x), where M[x) is the minimal 
polynomial of a over K and the denominator is m[ai,ai) = M'[ai). For each 
conjugacy class {0;^^ , . . . , ai - } of roots of M[x) over K(q;), we have that 

m[ai^,x) m[ai,,x) 

> — ? T^ifc [t) = trace— rUi^ [t). 

^^m[ai^,aij m{ai,,a^,) 

Where the trace is computed for the extension K[a,t,x) C K[a,t,x)[ai). 
Hence, we need to compute only one term of the Laurent interpolation for 
each conjugacy class of roots of M[x) over K[a). These conjugacy classes are 
determined by the factorization of M[x) in IK(a)[a;]. 
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Remark 10. To compute fast the trace of v — ^^^^^^u^ S K(a, Q!i)[x], first, 

we can compute the Newton sums 1 1 ' " ~ 1 from the minimal 

polynomial of ai-^ over K(a). If the coefficients of v are polynomials in t, we 
compute easily the trace of v computing the trace of each coefficient of v. If 
the coefficients of v are not polynomials in t, we can write v as n/{t + b), 
b G K(a,aj), n G K{a, ai)[t\. This is due to the fact that the variable t only 
appears on the term Ui and it is a linear fraction. Now, let g(t) be the minimal 
polynomial of —b over K(a) and gi(t) — g{t)/{t + b) E 'K{a,ai)[t]. Then v = 
njit + &) — (n ■ gi{t))/g{t) and trace{v) — trace{n ■ gi)/g(t) can be easily 
computed. 

Thus, we can compute the polynomial F (i.e. the standard parametrization) 
computing gcd and traces and norms in some fields of the form K(a, a^). To 
sum up, our algorithm to compute the standard parametrization of U is the 
following. 

Algorithm 11. Input: A curve C given by a proper parametrization ip{a,t) 
with coefficients in ]K(a) . 

Output: Either C is not defined over K or cj), the standard parametrization 
of the hypercircle associated to "tpit). 

1. Set M{x) the minimal polynomial of a over K. 

2. Set m{a,x) = M{x)l{x - a) e K{a)[x]. 

3. Compute m{a, x) — fi{x) ■ ■ ■ fr{x) the factorization ofm{a, x) over K(q:). 
I SetF=^teKia,t)[x]. 

5. For 1 < i < r do 

(a) Set OLi a root of fi(x). 

(b) Set = 'ijj{ai,t) the parametrization of the curve C"'^ . 

(c) Compute three good parameters ti,t2,t^ in the sense of remark^ 

(d) If two parameters ti, tj are found .such that ip{ti) and ^{tj) are well 
defined but not attained by 7/;'^' then Return C is not defined over K. 

(e) Compute Sk such that ip{tk) = ip''^{sk), 1 < fc < 3. 

(f) Compute Ui{t) = the linear fraction such that u{tk) = s^. 

(g) If tp ^ Tp'^^lui) the Return C is not defined over K. 

(h) Compute v — m{ai, x) /m{ai, ai) ■ Ui{t) G ¥^{a,t,ai)[x\. 

(i) Compute w = trace(v) for the extension K(a,i,x) C M^(a,t,x){ai). 
(j) Set F ^ F + w eK{a,t)[x]. 

6. Write F = (j}o{t) + 0i(t)a; + . . . + (?!)„_i(t)a;"-i. 

7. Return 4> = (00, • ■ • , 0n-i)- 
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Example 12. Now we present a full small example of the algorithm. Let K = Q, 
a a root of M{x) = — 2, consider the proper parametrization ip of a plane 
curve: 

X = {{lla^ + 15a'^+9a+ll)t^ + {7a^ + Ua'^ + l4a+7)t'^ + {a^ + 2a'^+4:a+l)t)/D 
y = ((15a3 + 9a2 + Ua + 22)t^ + {2ba^ + 29a^ + 16a + 25)t^ + {9a^ + 18a^ + 
15a + 9)t + a^ + 2a^ + 4a + 1) jD, 

withD = {7t^ + {12a^+3a'^+6a+l2)t'^+{6a^+l2a'^+3a+6)t+a^+2a'^+Aa+l). 

Now, M{x) = (x — a)(a;+a)(a;^+a^) is the factorization of M{x) in'K{a)[x\. 
m{a,x) = {x + a){x^ + a^) and m{a,a) = 4a^ = M'(a). start with F = 
^g^t = l/SiaxH + a^xH + a'^xt + 2t). 

From the factors of m{x) we have two conjugacy classes of roots of m over 
K(a). The first one is {—a}- Let a be a Q^- automorphism such that cr(a) = —a. 
Hence, we consider the conjugate parametrization ijj'^ : 

X = {(-lla^ + ISa^ -9a + ll)t^ + {-7a^ + Ua^ - Ua + 7)t^ + (-a^ + 2a^ - 

Aa + l)t)/Di, 

y = ((-15a3 + 9a^ - lla + 22)t^ + {-2ba^ + 290^ - 16a + 25)^^ + {-9a^ + 
18a2 - 15a + 9)t - a^ + 2a^ - 4a + l)/Di, with Di = {7t^ + (-12a3 + Sa^ - 
6a + 12)^2 + (-6a3 + 12a2 - 3a + Q)t - a^ + 2a^ - 4a + 1). 

We have to compute the automorphism such that tjj^t) = ij)'^ {u„{t)) . we 
evaluate {ip'')~^{ijj{tk)) and obtain: 

m = rio) 

= ^<^(8/31a^ - 4/31a2 + 2/31a - 1/31) 

V'(2) = V'^(128/511a^ - 32/511a^ + 8/511a - 2/511) 

Hence, = is such that u„{Q) = 0, u^{l) = 8/31a^ - 4/31a^ + 2/31a - 
1/31, Ua{2) = 128/511a3 - 32/511a'2 + 8/511a - 2/511. We can com,pute 
u(t) = by solving a linear homogeneous system of equations and get the 

solution 

, , a^t 

In this case tp'^{urj) = so C = . We can update F by adding: 

, N , / \ / \ —x^t + ax^t — a^xt + a^t 
m{-a,x)/m{-a,-a)ua{t) = ^ 

So now: 

_ ax^t^ + a'^x'^t^ + ax'^t + a^xt^ + 2t^ + a^t 
~ 8< + 2a3 

For this root, all operations are done in K(a) since cr(a) = —a e K(a). 

Now, we have to deal with the roots ofx^+a^. Let (3 be a root ofx^+a^ andr 
a Q- automorphism such that T{a) = (3. Consider the conjugate parametrization 
-tjjr: X = (((-lla2 + 9)/3-15a2 + ll)t3 + ((-7a2 + 14)/3-14a2 + 7)t2 + ((-a2 + 
4)/3-2a2 + l)t)/L>2, y = (((-15a2 + ll)/3-9a2 + 22)i3 + ((-25a2 + 16)/3-29a2 + 
25)t2 + ((-9a2 + 15)/3 - ISa^ + 9)t + (-a^ + 4)/3 - 2a2 + l)/D2, where D2 = 
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7t3 + ((-12a2+6)/3-3a2+12)i2 + ((-6a2+3)/3-12a2+6)t+(-a2+4);9-2a2 + l. 
In this case, we are taking the relative base [a^/S^ | < i < 3, < j < 1} o/ 
Q(a,/3) over Q. Now we compute Ut such that 4'it) = ip'^{u{t)). for this 

V'(0)=^^(0) 

= V^^((2/9a2 - 2/9a + 1/9)^ + 2/9a^ - l/9a + 1/9) 
^{2) = (32/ 129^2 - 16/129a + 4/129)/3 + 32/129a^ - 4/129a + 2/129) 
From this data, we can compute: 

t 

Ur{t) 



{a -/?)< + 1 



//7 is the other root of +Q:^ ( i.e. j — —j3) and 5 is a 'Q^- automorphism such 
that 5{a) = 7, then Uy{t) — t/{a — "/)t + 1. We have to compute the trace of 

m(l3, x) , , Pix^t - o?x^t - cx^fixt + 2t 

V — U it] 

m(/3,/3) ' (-8/3 + 8a)t + 8 

over Q{a,t). This is done using the technique described in Remark \l(A 
'O^x^t^ - a^x^t^ - a^x^t + 2xt^ + 2at^ + 2t 



trace{v) 



8aH^ + Sat 



To compute this part, we have made computation in Q{a,f3). We add trace{v) 
to F and get 



where 



F — (j)iX + (j)2X^ + 4>3X^ 



2t^ + 3a^t^ + ia^t^ + at a^t^ + 2a^t^ + at^ 



^ 8^3 + 6aH^ + AaH + a ' ^ + 6aH^ + AaH + a ' 

_ a^f^ + at^ _ af^ 

^'^ ~ %t^ + &aH^ + 4,aH + a''^^ ~ 8t^ + 6aH^ + AaH + a 
And (j) = {(f>o, 01, 02, 03) the .standard parametrization of the hypercircle asso- 
ciated to -0. 

So far, Algorithm [TT] only computes the hypercircle U. The algorithm is 
able to detect if C is not defined over IK, but apart from that it does not provide 
much more useful information. In the rest of the section, we show that, if C 
is not defined over K, how can we compute the minimum field L such that 
K C L C K{a) and C is defined over L. Note that K(a) always is a field of 
definition of C, so the existence of L is always guaranteed. 

Theorem 13. Let C be a curve not K-definable but 'K{a)-parametrizable. Let 
L be the minimum field of definition of C containing K. K C L C K(a). Then 
L is the subfield of the normal closure K(a) over K that is fixed by the K- 
automorphisms a o/K(a) such that C ~ . 
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Proof. First, we recall that the intersection of fields of definition of C is a field 
of definition of C. Hence, since IfC(a) is a field of definition, there always exists 
a minimum field of definition L of C containing K. 

From [2] it follows that if Li C L2 is any algebraic finite normal extension 
and L2 is a field of definition of C, then Li is a field of definition of C if and 
only if C = for all a e ^^^(L a/Li). 

L et G ^ {a e Aut(K{a)/K) \ C = C}. Clearly, G is a subgroup of 
Aut{K{a)/ K). T his follows from the f act th at {Cy = 0'"°" . Let L be the 
subfield of K(a) that is fixed by G. L C K(a) is a normal extension and, if a is 
a L-automorphism of K(a) then ct € G so C = C^. In this conditions, L is a field 
of definition of C. Moreover, it is the smallest field of definition of C containing 
K. If IK C Li C L is a subfield of L, then Gi, the set of Li-automorphisms 
of K{a), is Gi 2 G. Hence, there is an automorphism t € Gi \ G. But then 
C ^ and Li cannot be a field of definition of C. Now, since K(a) is also a 
field of definition of C, then L C K(q;). □ 

If (To = Id,ai, . . . ,an-i are the automorphisms defined in Section [21 then 
for any a G Aut{K{a)/K), it happens that = C^* for some i, 0<i<?i — 1. 
Hence 

L= Pi {x e K{a) \ ai{x) = x} 

0<i<n-l 
0=0"^ 

If C is not defined over K, we compute in step 5 of Algorithm [Tl] the set of 
automorphism ai such that C — C"^ . For any such i, let m be the degree of 
over K(q;). If a; € IK(a), we can write (7i(x) = X]j=o ~ ^i'^h where U € K(a). 
X is CTi invariant if and only ii x — Iq, li — 0, 1 < z < m — 1. This provide a set 
of K-linear equations in the coordinates of x in K{a) = K". Note also that if 
tti and ttj are conjugate over IK(q;), the equations imposed by ai and aj are the 
same. Hence, we only need to compute them once for each set of conjugate roots 
of M{x) over K(a). Solving the system of linear equations provide a base of L 
as a K-subspace of K(q;). From this equation, we may reapply Algorithm 1111 
but to the extension L C K(q;). In this case we already have computed the 
automorphisms so we can reuse this computation. 

4 Complexity and Running Time 

We now compute the complexity of Algorithm [11] in terms of number of op- 
erations over the ground field K. The analysis is by no means sharp, we only 
intend to prove that there is a polynomial bound and that the main obstacle is 
the degree of a over K. 

Theorem 14. Let K fee a computable field with factorization of characteristic 
zero, a algebraic of degree n over K of minimal polynomial M{x). Let 'ip{t) = 
(V'Oj ■ • ■ : V'm-i) be a proper parametrization of a spatial curve C with coefficients 
in K(a). Then the number of operations over K of Alaorithm \ll\ is bounded by 
K + 0{md^n^) where K is the time needed to factor M{x) in K{a)[x]. 
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Proof. We only use naive algorithms. The factorization of M[x] can be per- 
formed standard methods [4l[T3| from a factorization algorithm in K[x]. Addi- 
tion in K(a) costs n operations and multiplication costs 0(n'^) operations and 
inversion O(n^). If /? is a conjugate of a, the worst case complexity of addition 
in K(a,/3) is O(n^) while multiplication is 0{n'^) and inversion O(n^). If / and 
g are two polynomials of degree at most d, their gcd costs 0{(Pn'^ +n^d'^) oper- 
ations in K{a) or 0{d^n'^ + n^d^) if their coefficients live in K(a,/3). Steps 1 — 3 
of the algorithm cost K + 0{n^). Step 4 is evaluating a polynomial in K(a), 
invert the result and multiply the polynomial this result. By Horner's method it 
is 0{n^). Step 5.6 can be done in 0{dmn) operations. For a parameter tk doing 
steps 5.C — d is evaluating m rational functions in K(a) and then compute m — 1 
gcd in K(a,/3) of degree d, this costs 0{md^n^ + mn^d'^). From Thcorcm|6]we 
have to try at most 0{d^ +n) times, so the total cost is bounded by 0{md^rJ). 

Computing step 5./ is just solving a system of 3 linear equations in 4 un- 
knowns in K(a,/3). This can be done in 0{n'^) operations. Now, comparing tj} 
and ij)" {u) in 5.e can be done evaluating both functions in 0{d) parameters. 
Each evaluation costs 0{n^), so in total, this step can be done in 0{mdn^). 
Step b.h we already have m(ai,x) lm{ai,ai) precomputed by conjugation, so 
we only need to multiply the polynomials, which is dominated by computing 
Oin) products [u is always of degree < 1 and we do not need to do anything 
with the denominator). This costs 0(n^). Now, instead of computing the min- 
imal polynomial of the pole of m, we can compute its characteristic polynomial 
over K(a). Since the characteristic polynomial of an n x n matrix can be done in 
0{rL^) operations and the matrix will have entries in K(a), we can compute this 
characteristic polynomial in 0{n^) operations. Step 5.j can be done in O(n^) 
operations. Hence step 5 is bounded by 0{md^'n7). Since we have to perform 
step 5 at most n times. We get a bound of 0{md^n^) operations over IK. 

If C is not K-definable. In step 5 we compute the automorphisms ai such 
that C — . From this automorphisms, we can compute the field of definition 
L in 0{n'^) operations and repeat the whole algorithm. It is clear that the 
running time for the extension L C K(a) is bounded by the case IK C K(a). So 
the global bound does not change. □ 

This result agrees with experimentation, the most important parameter is 
the degree of a over IK and the ambient dimension of C tend to be not relevant 
in the algorithm compared to the other parameters. 

Computing the hypercircle using Definition [1] is too slow, because we have 
to work with an ideal in n variables over K and make the quotient by the ideal 
defined by the denominator. In [8] the authors proposed a method to compute 
the parametrization of the hypercircle. It is based in the following result. 

Theorem 15. Let ^ be a proper parametrization of C with coefficients in 
K(q;). Let G{xi, . . . ,Xm) : C ¥ be the inverse of the parametrization. G € 
K{a){xi, . . . , Xm). Write G = X)r=o ^i^^ ' ^ ^{xi, . . . , x^), < i < n — 1. 
Consider (p = (Goiip), . . ■ , G„_i('(/')). Then C is defined over K if and only if <f> 
is well defined and parametrizes a curve in F". In this case cj) is the standard 
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parametrization of the associated hypercircle to ip. 
Proof. See [5] 



□ 



Algorithni[TT]and the algorithm in Theorem llSl have been implemented in the 
Sage CAS [TU] , the code for the method presented in this paper can be obtained 
from [TT]. We are interested in the average case, so we will assume that our 
curve is planar (since we can always make a generic projection). However, the 
method presented here also performs well for spatial curves. For the method 
based on the inverse of the parametrization of [8] we compute (Go, . . . , G„_i) 
but we do not simplify the composition Gi^ip). This is done to avoid artifacts 
in the running time that appeared if we simplify the composition. The inverse 
of ip is computed using the resultant method explained in [5] . 

We show the results for random curves of degree 2, 5, 10, 25 and 50. First 
over an extension of Q of degree 2, over a cyclotomic extension of degree 6 and 
a random extension of degree 5. In all these cases C is defined over K. 



Case: + 1 = 


method \ degree of C 


2 


5 


10 


25 


50 


Moving hyperplanes; 
Inverse-based method: 


0.08 
0.03 


0.15 
0.15 


0.27 
13.16 


0.87 

> 60 


2.58 



Case: a'' + + a* + + + a + 1 = 


Moving hyperplanes: 
Inverse-based method: 


1.12 
0.13 


2.00 
28.61 


3.71 

> 60 


13.13 


48.01 



Case: a of degree 5, random minimal polynomial 


Moving hyperplanes: 
Inverse-based method: 


1.04 
0.12 


2.01 
10.36 


3.84 

> 60 


14.28 


> 60 



Now, we show a table with a random extension of degree 5 but C is not 
defined over K. In this case is more evident that the moving hyperplanes method 
is better. With high probability it will detect that the curve is not defined over 
K while trying to compute the automorphisms u(t), on the other hand, the 
inverse-based method always has to compute the inverse of the parametrization 
ip. In all cases, our algorithm computed the minimum field of definition of the 
corresponding curve. 



Case: a of degree 5, C not defined over K 


method \ degree of C 


2 


5 


10 


25 


50 


Moving hyperplanes: 
Inverse-based method: 


0.36 
0.08 


0.59 
12.68 


0.91 

> 60 


2.16 


5.13 
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